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ABSTRACT 

We have used multi-band high resolution HST WFPC2 and ACS observations 
combined with wide field ground-based observations to study the blue straggler 
star (BSS) population in the galactic globular cluster NGC 6388. As in several 
other clusters we have studied, the BSS distribution is found to be bimodal: 
highly peaked in the cluster center, rapidly decreasing at intermediate radii, and 
rising again at larger radii. In other clusters the sparsely populated intermediate- 
radius region (or "zone of avoidance" ) corresponds well to that part of the cluster 
where dynamical friction would have caused the more massive BSS or their binary 
progenitors to settle to the cluster center. Instead, in NGC 6388, BSS still popu- 
late a region that should have been cleaned out by dynamical friction effects, thus 
suggesting that dynamical friction is somehow less efficient than expected. As by- 
product of these observations, the peculiar morphology of the horizontal branch 
(HB) is also confirmed. In particular, within the (very extended) blue portion of 
the HB we are able to clearly characterize three sub-populations: ordinary blue 
HB stars, extreme HB stars, and blue hook stars. Each of these populations has 
a radial distribution which is indistinguishable from normal cluster stars. 



Subject headings: Globular clusters: individual (NGC6388); stars: evolution - 
Horizontal Branch - binaries: general - blue stragglers 
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INTRODUCTION 



The present work is a continuation of a series of papers exploring the interplay of 
dynamics and stellar evolution in Galactic globular clusters (GCs), through the detailed 
study of their peculiar stellar populations. To date we have primarily focused on Blue 
Straggler Stars (BSS). BSS are brighter and hotter than the Main Sequence (MS) Turnoff 
(TO) and are thought to form by the evolution of primordial binaries and/or by the effect 
of stellar dynamical interactions. They are most abundant in the central regions of GGCs, 
which are usually very congested in the optical bands, but relatively open in the UV. To 
fully understand the role of dynamics one must track how the population of BSS (or any 
other population arising from binary evolution or stellar interactions) varies over the entire 
radial extent of the cluster. The results of the earlier papers in the series are summarized in 
Table ffl 

Of the clusters we have explored so far NGC 6388 is perhaps the strangest. As shown 
for the first time by Rich et al. (1997), this bulge cluster has an extended blue H orizontal- 
Branch (HB), despite its high metal content ([Fe/H] = — 0.44; ICarretta et al.ll2007l ). Such an 
HB morphology is contrary to expectations: metal-poor clusters usually have predominantly 
blue HBs, while metal- rich clusters have stubby red HBs. NGC 6388 is a vivid example of the 



so-ca lled second-parameter phenomenon (see, e.g.. ISweigart fc Catelanlll998l ; ICatelan et al. 



20061 ). Moreover, hints for the presence of an intemediate mass black hole (IMBH) at the 
center of NGC 6388 have been recently inferred from the observed surface brightness profile 
(which shows a devi ation from a flat core behavior in th e inner ~ 1") and from detailed 



dynamical modeling (jMiocchi 1 120071 ; lLanzoni et al. Il2007d ). More robust evidence (through 



kinematics, X-ray and/or radio observations) is needed to confirm this finding. 

In this paper we present multi-wavelength observations of NGC 6388. We have combined 
high-resolution ultraviolet (UV) and optical images of the cluster center obtained with HST, 
with complementary wide-field observations from the ground, covering the entire cluster 
extension. The data and photometric reductions are described in § [5J A general overview of 
the color-magnitude diagram (CMD) is discussed in § [31 The BSS population is described 
in § HI and the Discussion is presented in § [51 



1 Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on WFI observations 
collected at the European Southern Observatory, La Silla, Chile. 
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2. OBSERVATIONS AND DATA ANALYSIS 
2.1. The data sets 

We have used a combination of high-resolution and wide-field photometric data sets: 

1. The High resolution set consists of a series of public multiband (from the UV, to 
the optical) WFPC2 and ACS HST images, which have been retrieved from the ESO/ST- 
ECF Science Archive. The WFPC2 images were obtained through filters F2b5W (mid-UV) 
and F336W (U), with total exposure times (t exp ) of 9200 and 1060 s, respectively (Prop. 
8718, P.I. Piotto), and trough filters FA39W (B) and FhhhW (V), with t exp = 370 and 62 
s, respectively (Prop. 6095, P.I. Djorgovski). This combined dataset allows us to examine 
both the hot (HB and BSS) and the cool (red giant branch and subgiant branch, RGB and 
SGB, hereafter) stellar populations in the cluster. In this dataset the planetary camera (PC, 
with the highest resolution of ~ 0'.'046 pixel -1 ) was roughly centered on the cluster center, 
while the wide field cameras (WFCs, at a lower resolution of ~ O'.'l pixel -1 ) sampled the 
surrounding outer regions. The photometric reduction of these images was performed using 
ROMAFOT (Buonanno at al. 1983), a package developed to perform accurate photometry in 
crowded regions and specifically optimized to handle under-sampled Point Spread Functions 
(PSFs; Buonanno & Iannicola 1989), as in the case of the WF chips. The ACS dataset is 
composed of a series of images (Prop. 9821, P.I. Pritzl) obtained through filters F435W 
(B) and F606W (V), with t exp = 11 and 7s, respectively. It gives complete coverage of 
the central cluster region out to 110" from the center (see the following section). All the 
ACS images were corrected for geometric distortions and effective flux (Sirianni et al. 2005). 
The photometric analysis was performed independently in the two images by using the the 
aperture photometry code SExtractor (Bertin & Arnouts 1996) and adopting an aperture 
radius of 3 pixels (corresponding to 0'.'15). 

2. The Wide field set is a complementary set of public B and V images obtained with 
the Wide Field Imager (WFI) mounted at the 2.2m ESO-MPI telescope at La Silla ESO- 
Observatory and retrieved from the ESO Science Archive. The WFI is a mosaic of 8-CCD 
chips, each of 2000 x 4000 pixels, with a pixel size of ~ 0'/24. It has exceptional image 
capability, with a global field of view (FoV) of 33' x 34'. Thanks to such a FoV, this dataset 
covers the entire cluster extension with the cluster roughly centered on the CCD #7. The 
WFI images have been corrected for flatness, bias and overscan using I RAF tools. Th e PSF 



fitting was performed independently on each image using DAOPHOT (jStetson! 119871 ). 
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2.2. Astrometry, Photometric Calibration, and Sample Definition 

Using the procedure described in Ferraro et al. (2001, 2003) the WFI catalogue has been 
placed on the absolute astrometric system. The 8 WFI CCDs have been astrometrized by 
cross-correlating each of them with the new astrometric 2-MASS catalogue using a specific 
tool developed at Bologna Observatory. Several hundred astrometric reference stars were 
found in each WFI chip, thus allowing an accurate absolute positioning of the sources. As a 
second step, a few hundred stars in the overlapping area between the WFI, and the WFPC2 
and ACS FoVs were used as secondary astrometric standards, in order to place the HST 
catalogs on the absolute astrometric system. At the end of the procedure the rms residuals 
(that we take as representative of the astrometric accuracy) were of the order of ~ 0'/3 both 
in RA and Dec. 

By using the procedure described in Ferraro et al. (1997, 2001), the photometric cal- 
ibration of the UV m agnitudes (F2bbW and F336W) has been performed in the STMAG 



system, adopting the iHoltzman et al. I (119951 ) zero-points. The optical (B and V) magni- 



tudes hav e been transformed t o the Johnson system by using the stars in common with the 



catalog of iPiotto et al. I (120021 ) . Linear transformations were adopted, and only small color 



equation terms were required to correct the response of the different filter profiles. 

Final lists with the absolute coordinates and homogeneous magnitudes for all the stars 
in the three considered catalogs were obtained. To minimize incompleteness effects in the 
ground based observations of the crowded central regions of the cluster, while still taking 
advantage of the superior capability of UV observations in detecting the BSS (Ferraro et al 
1999a, 2001), we divided the dataset in two main samples: The HST sample, includes only 
stars in the WFPC2 and complementary ACS catalogs. It covers approximately the inner 
r < 110" of the system (with the WFPC2 FoV almost entirely included in the ACS FoV; see 
Figured]). The WFI sample includes only stars observed with WFI and lying beyond the 
WFPC2 and ACS FoVs. The WFI sample covers the outer regions, well beyond the cluster 
extension (see Figure [2]). 

By combining the data sets described above, with additional images of the cluster center 
obtained with the HST ACS High Resolution Camera, very accurate determinations of the 
center of gr avity, surface density p rofile, and surface brightness profile have been recently 



P 

7c) 



obtained by lLanzoni et al. I (l2007d ). In particular, it has been found that the observed pro- 



files show a deviation from a flat core behavior in the inner ~ 1", suggesting that NGC 6388 
might host an IMBH of ~ 5.7 x 1O 3 M in its center. However, by excluding the points at 
r < 1", the density profile is well fit by an isotropic single-mass King model. The result- 
ing cluster structural parameters (concentration, core radius and tidal radius) are listed in 
Table [2] together with the new estimate of the center of gravity. 
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3. CMD OVERVIEW 

3.1. The HST Sample 

The CMDs of the HST sample in the UV and the optical bands are shown in Figures [3] 
andHl respectively. As apparent, all of the cluster evolutionary sequences are clearly defined 
and well populated. 

Particularly notable is the Horizontal Branch (HB) morphology. Beside the red clump, 
which is a typical feature of metal rich stellar populations, the HB clearly shows an extended 
blue tail (BT), first noticed by Rich et al. (1997) and by Piotto at al. (1997). Among a 
total of 1763 HB stars counted in the HST sample, five sub-populations can be distinguished 
(see Sect. 13.1.11 and 13.1.21 for details): (i) the red-HB (RHB) population, consisting of 1418 
stars grouped in the red clump; (ii) 15 RR Lyrae variables, which we identified by cross 
correlating the positions in our catalog with those published by Pritzl et al. (2002) □; (Hi) 
267 blue-HB (BHB) stars ; (iv) 26 Extreme-HB (EHB) stars ; and (v) 37 Blue Hook (BHk, 
to avoid confusion with BH for black hole) stars. 

Several previous works have shown that the HB morphology in NGC 6388 is complex. 
A new extensive study, based on much of the same observational data used h ere, and also 



discus sing the HB morphology of NGC 6441, has recently been published by iBusso et al. 



( 120071 ). In the present paper we take advantage of the complex HB structure of NGC 6388 
to review the HB nomenclature, which has become rather confused in the literature and is 
often ambiguously used. Then, we briefly discuss the blue HB sub-populations and the HB 
red clump, the latter being used as cluster reference population for the study of the BSS 
radial distribution (Sect. I4T2]) . 



3.1.1. The HB population: nomenclature and radial distribution 

The HB is composed by helium-core/hydrogen-shell burning stars. It is traditionally 
split into red, variable, and blue (RHB, VHB, and BHB, respectively), depending on whether 
the stars are redder than, within, or bluer than, the RR Lyrae instability strip. 

The concept of HB blue tails probably originated with the CMD of NGC 6752, which 



2 Since our photometry is just a snapshot, the position of each variable star in the CMD is not an 
indication of the mean properties. The remaining stars found within the "RR Lyrae region" of the CMD, 
but not included in the Pritzl et al. catalog are not considered in the following analysis, since they possibly 
are field contaminating stars. 
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was presented by Russell Cannon at the 1973 Frascati globular cluster workshop, but not 
published for many years. In visual CMDs of NGC 6752 (and many others to follow), the 
HB drops downward at high temperature often becoming an almost vertical sequence. This 
feature looked li k e a ta il hanging from the hor izonta l part of the BHB, hence the name. 
Rood fc Crocker! (Il989l ) , iFusi Pecci et al.l (119931 ) , and iRecio-Blanco et al.l (120061 ) have each 
suggested ways to measure BTs, and the fact that measures of BTness keep being invented 
demonstrates a lack of consensus on a definition of BTs. In addition, sub-populations like 
EHB or BHk stars are sometimes recognized within the observed BTs, even if without a 
precise observational definition. 

The extreme HB population is theoretically well defined: EHB stars lie at the hottest 
extreme of the zero-age HB (ZAHB), and they do not return to the asymptotic giant branch 
(AGB), but rather spen d their He-shell bur ning phase as hot AGB- manque or Post-early 
AGB stars (for example, iDorman et al.lll993l ). There is no comparably precise way to obser- 
vationally select EHB stars. If far-UV (FUV) (e.g., HST F160W) photometry is available, 
detailed comparisons with stellar models can be made. These suggest that in a few clusters, 
the trans ition between BHB and EHB stars may be associated with a gap in the HB mor- 
phology (IFerraro et al.lll9981 ). In the present paper we have assumed this to be the case (see 
below), but at this point that is only an assumption. The importance of EHB stars is also 
connected with the fact that they and their progeny are thought to be the source of the UV 
radiation excess observed in the integrated spectra of some elliptical galaxies (Dorman et al 
1995; Yi et al. 1998), and one might be able, for example, to determine the age of the galaxy 
on the basis of its UV excess. In this context NGC 6388 plays a particularly important role, 
since it is one of the most metal-rich systems containing EHB stars that can be individually 
observed. 

In a few clusters, including NGC 6388, there is an additional population hotter and 
less luminous than the EHB stars. Following nomenclature used in recent studies, we call 
this population Blue Hook stars. In visual and even some UV (e.g., 717,255, — rn^o) 

CMDs, BHk stars appear as fainter extension of the BT and are separated from the EHB 
population by a gap. While the effective temperature T eff of HB stars can be reasonably well 
determined from their position along the BT, it is not appropriate to extrapolate this to the 
BHk. Accurate stellar parameters for BHk stars require FUV ph otometry (see for example 
the BHk studies in NGC 2808 and u Cen bv lMoehler et al. Il2004h . Indeed, it is only in FUV 
CMDs that the origin of the name "blue hook" becomes apparent. 

2007bl). 



Not all BT clusters have EHB stars (see the case of NGC 1904 Lanzoni et al. 



and not all clusters with EHB stars have BHk stars (see the cases of M13 and M80 



Ferraro et al 



19981 ). In order to clearly show the difference between cluster with BTs populated up to the 
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EHB region, and clusters with BHk stars, UV photometry is essential. In Figure [5] we com- 
pare the (m 255 , m 255 - U) CMDs of NGC 6388 to that of M80 (Ferraro et al. 1997, 1999a). 
We chose M80 because, among the clusters for which we have a full range of data, it is the 
one with the HB extending to highest T e g. From comparison with evolutionary tracks in 
the (mi6o, ^160 — V) plane (Dorman and Rood, unpublished) we know that the M80 HB is 
populated all the way to the extreme blue-end of the Z AHB. In this plane the HB of M80 
shows a clear gap at the transition from BHB to EHB (IFerraro et al.lll998l ). In Fig. [5] this 
gap is also visible at m 255 = 19. Since there is a corresponding gap in NGC 6388, we ten- 
tatively identify the stars with 19.0 < m 255 < 19.8 as EHB. The HB sequence of NGC 6388 
is significantly more extended than that observed in M80, where there is no analogous pop- 
ula tion at 777.255 — 20. For that reason we identify the latter as BHk stars in NGC 6388, as 



do iBusso et al.l (120071 ). The BHk extends well below the cluster TO in the optical CMD: 
V > 21 in the (V, B — V) plane (see Fig. Ej). Because of the uncertainties in detection 
of these faint stars using optical wavelengths in such crowded regions, our BHk sample is 
selected from the UV CMD (WFPC2 sample). By using these criteria we have defined the 
selection boxes sketched in Figs. [3] and El and obtained the values quoted in the previous 
section for the number of HB stars belonging to each sub-population. 

It has been suggested that EHB might originate in binary systems (Bailyn 1995), or 
they were formed through a collisional channel. Indeed, similar stars in the field and in 



old open clusters have been found to often belong to binary systems (I Green et al.l 12001 



Heber et al.l 120021 ) . In contrast, a recent study by Moni Bidin et al. (2006) has found no 
evidence of binarity in the EHB population of two globular clusters (M80 and NGC5986). 
The nature of BHk stars is still a m atter of debate. They may be r elated to the so-ca lled late 
hot flashers (IMoehler et al. 112004 ) . or to high helium abundances ( IBusso et al.ll2007l ). Given 
that the origin of EHB and BHk stars is still uncertain, it is useful to check whether their 
radial distributions show any suggestion of binarity or stellar interaction, as it is the case for 
BSS. We have therefore compared the cumulative radial distributions of the BHB, EHB, and 
BHk stars to that of RHB stars, which are representative of normal cluster populations (see 
Sect. I4.2p . As shown in Figured all radial distributions are consistent with being extracted 
from the same parent population, with Kolmogorov-Smirnov (KS) test probabilities of 59%, 
46% and 60%, respectively. Thus, the radial distribution of BHB, EHB and BHk populations 
is consistent with that of normal cluster stars, in agreement with previous findings by Rich 
et al. (1997), and possibly suggesting a non-binary nature for these systems. However, a 
binary EHB star could consist of a 0.5 M He-burning star with a 0.2 M & He white dwarf 
companion, i.e., a total mass smaller than the TO mass and comparable to RHB masses. As 
a consequence, even if the initial binary mass were large enoug h to have sunk to the cluster 
core, the central relaxation time of NGC 6388 is 8.3 x 10 7 yr (jDjorgovski Ill993l ). less than 
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the HB lifetime, so that an EHB binary could move outward, instead of being segregated 
into the center. 



3.1.2. The red HB clump and the distance of NGC 6388 

Since the HB red clump in this cluster is very well defined in the optical CMDs, we 
have selected the RHB population in this plane, and then used the stars in common between 
the optical and the UV WFPC2 samples to identify it in the (777.255, 777255 — U) CMD. The 
selected stars are marked as pentagons in Figs. [3] and El 

With our high-quality data set and such a well defined HB red clump, we have estimated 
the distance modulus and the reddening of NGC 6388 by comparing its CMD and luminosity 
function to those of the proto-type of metal-rich GCs: 47 Tuc. As shown in Figure [HJ other 
than the blue HB, the overall CMD properties of NGC 6388 closely resemble those of a normal 
metal-rich cluster. In order to overlap the CMDs, and align the HB red clump and the RGB 
bump of the two clusters, the color and the magnitude of NGC 6388 have to be shifted by 
5(B - V) = -0.28 and 5V = -3.15, respectively. Thus, by adopting (m - M) = 13.32 
and E(B - V) = 0.04 for 47 Tuc (Ferraro et al. 1999b), we obtain (777 - M) v = 16.59 
and E(B — V) = 0.32 for NGC 6388. This yields a true (unreddened) distance modulus of 
(777 — M)o = 15.60 ± 0.2 , which corre sponds to a distance of 13.2 ± 1.2 Kpc (to be compared 



with 10 Kpc quoted by lHarris Ill996l ). 



The proper comparison of the CMDs of the two clusters deserves additi onal comments- 



First , NGC 6388 is known to be sligh tly more metal-rich (|Fe/H| = —0.44; ICarretta et al. 



20071 ) than 47 Tuc ([Fe/H] ~ -0.6; ICarretta et all 120041 ) . However, current theoretical 



models suggest that such a small overabundance (5[Fe/H] ~ 0.15 dex) would generate just 
a small difference in the HB absolute magnitude (5My B ~ 0.03 mag). Second, the presence 
of differential reddening of the order of 0.07 (see Busso et al. 2007) can spread the HB 
red clump, increasing the uncertainties of the entire procedure. However, this contribution 
is significantly smaller than the conservative estimate of the error in the derived cluster 
distance, 8{m — M) = ±0.2 mag. Finally, the morphology of the HB red c l ump is not 



exactly the same in the two clusters. However, as discussed by ICatelan et al.l (120061 ). who 
compared the two CMDs by using a reddening-independent quantity, the main difference 
between the two red HB clumps consists in the fact that the bluer RHB stars in NGC 6388 
(~ 20% of the total RHB population) are slightly brighter than the average in 47 Tuc. This 
feature might be interpreted in the framework of a sub-population with a higher helium 
content. However comparison of the luminosity functions in Figure [8] clearly shows that the 
relative position of the HB red clump and RGB bump is quite similar in the two clusters. 
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Since the location in magnitude of the RGB bump is quite sensitive to the helium content 
(see Fusi Pecci et al. 1990), such a nice correspondence clearly demonstrates that at least 
the main component of the stellar population of NGC 6388 has an helium abundance fully 
compatible with that of 47 Tuc, while only a minor fraction of the cluster stars could be 
helium enhanced (this is also in agreement with the findings of Catelan et al. 2006, and 
Busso et al. 2007). The possible impact on the relative distance of the two clusters derived 
above is therefore negligible. 



3.2. The WFI Sample and Background Contamination 

Figure [9] shows (V, B — V) CMDs of four radial zones of the WFI sample. The se- 
quences seen in Fig. H] are still obvious in the interval 120" < r < 250", but there is also 
significant contamination from field stars. The CMDs are progressively more contaminated 
as r increases: cluster sequences are barely visible in the region 250" < r < 490", less so 
for 490" < r < 800", and have vanished for r > 800". The contamination has two main 
components (the bulge and the disk populations of the Galaxy): the first is an almost ver- 
tical blue sequence with 0.5 < (B — V) < 1.0, the second is another vertical sequence with 
(B — V) ~ 1.3, which clearly indicates the presence of metal rich stars. 

Indeed, Figure M shows that field contamination is particularly severe in this cluster. 
For this reason we decided to limit the following analysis to the HST sample, and to use 
the most external region of the WFI sample (r > 800") to statistically estimate the field 
contamination level. We have counted the number of background stars and derived an 
appropriate background density for each selection box discussed in the paper (see Sect. @] 
for the definition of the BSS population). Then, we have used this background density to 
statistically decontaminate e ach population in the HST samp l e, by following a procedure 



similar to that described by iBellazzini. Ferraro &: Buonanno I (119991 ) . Table H] shows the 



observed sample size and the resulting statistical estimate of field contamination for each of 
the sub-populations in each radial region of the cluster. Statistical decontamination has the 
disadvantage that we don't know whether a given star is a member of a given sub-population 
or not. However, all of the conclusions of the current paper depend on number counts, so 
that background correction only increases the noise without affecting the conclusions. Future 
proper motion studies currently ongoing for this cluster will finally assess the real membership 
of each star. Preliminary results indicate contamination counts which are in agreement with 
those listed in Table 4. 



-10- 



4. THE BLUE STRAGGLER STAR POPULATION 

4.1. The BSS Selection 

Hot populations like BSS and BHB stars are the brightest objects in UV CMDs, while 
the RGB stars, that dominate the emission of GCs in the optical bands, are faint at these 
wavelengths. In addition, the high spatial resolution of HST minimizes problems associated 
with photometric blends and crowding in the high density central regions. Thus HST UV 
CMDs are the optimal tool for selecting BSS in GCs. Given this, our main criterion for 
the selection of the BSS population is the position of stars in the (m 2 55, "^255 — U) CMD. 
To avoid incompleteness biases and contamination from TO and sub-giant stars, we adopt a 
magnitude threshold that is about one magnitude brighter than the TO point: m 2 55 = 21.85. 
Figure [3] shows the adopted selection box and the candidate BSS in the UV CMD. Using 
the BSS in common between the WFPC2 and the ACS FoVs, we have transformed the BSS 
selection box from the UV plane into the optical plane. To avoid regions with very high risk 
of Galactic contamination, we have considered only stars with (B — V) < 0.7. The resulting 
candidate BSS in the complementary ACS field are shown in Fig. [61 

The final sample is of 153 BSS in the HST sample: 114 are found in the WFPC2 dataset, 
and 39 in the complementary ACS sample. The magnitude and the positions of the selected 
BSS are listed in Table [31 which is available in full size in electronic form. 

4.2. The BSS Projected Radial Distribution 

In order to study the radial distribution of BSS (or any other population) for detecting 
possible peculiarities, a reference population representative of normal cluster stars must be 
defined. 

In our previous papers we have used both the RGB or the HB as reference populations. 
In NGC 6388 the RGB population is affected by a significantly larger field contamination, 
with respect to the HB. On the other hand, the HB morphology is quite complex, the presence 
of a BT in such a metal-rich cluster is unusual, and the nature of EHB and BHk stars is 
still unclear (see Sect. 13.1.11) . Instead, the HB red clump is a common feature of similar 
metallicity GCs, it is bright and well defined both in the UV and in the optical CMDs, and 
it comprises the majority (80%) of the HB population (see § 13.11) . We have verified that RHB 
and RGB stars share the same radial distribution over the region (r < 110"), suggesting that 
RHB stars are indeed representative of normal cluster populations. For all these reasons, we 
have chosen the RHB as the reference population for NGC 6388. 
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We first compare the BSS and the RHB cumulative radial distributions. Since we expect 
a negligible number of field stars contaminating the BSS and RHB sample (see Table 4) no 
correction has been applied to the observed sample used to construct these cumulative radial 
distributions. As shown in Figure [TDl the trend is bimodal, with the BSS more segregated 
than RHB stars in the central cluster regions, and less concentrated outward. This result 
replicates what we fo und in a numbe r of other clusters a nd listed in Table 1 (See also for M5 



( jWarren et al. 1120061 ) and for M55 (IZaggia et al.lll997l )). The KS probability that the two 



populations are extracted from the same parent distribution is ~ 10 4 . 

For a more quantitative analysis we have divided the surveyed area in 6 concentric 
annuli (sketched in Fig. [Tj), and the number of BSS (iVess) and RHB (AT RH b) stars was 
counted in each annulus (see the values listed in Table HJ). We have then computed the 
double normalized ratio (Ferraro et al. 1993): 



R 



N /N tot 



pop 



P°P £samp /L S ^ P ' 

where pop=BSS, RHB, iVp 0p refers to statistically decontaminated number counts (see 
Sect. I3.2p . and the luminosity in each annulus has been calculated by integr ating the single- 



mass King model that best fits the observed surface density profile (see lLanzoni et al 



2007d ). with the distance modulus and the reddening previously quoted, and by properly 



taking into account the incomplete spatial coverage of the outermost annulus. 

As expected from stellar evolution theory (Renzini & Fusi Pecci 1988), the radial trend 
of -Rrhb is essentially constant with a value close to unity. On the contrary, the BSS radial 
distribution is very different and is clearly bimodal. As shown in Figure [TTJ -Rbss reaches 
a value of almost two at the center, while -Rrhb has no central peak. -Rbss decreases to a 
minimun near r = 5r c , and rises again near r = 11 r c . 



5. DISCUSSION 



In Figure [121 we have compared the radial distribution of the ratio between the BSS and 
HB number counts computed for N GC 6388, with that obtained for other GCs showing a 
bimodal distribution (see for example lLanzoni et al.ll2007l ). The position of the observed min- 
imum in NGC 6388 resembles that of M3, but NGC 6388 has a core radius ~ 3 times smaller. 
In physical units its minimum is closer to the cluster center than in any previously observed 
cluster. By equating the dynamical friction timescale (tdf oc a 3 / p\ iBinney &: Tremaine 1119871 . 
see also Mapelli et al. 2006) to the cluster age (assumed to be t = 12 Gyr), one can estimate 
the value of the radius of avoidance (r avoi( j). This is defined as the radius within which all 
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the stars of ~ 1.2 M Q (the expected average mass for BSS) have already sunk to the core 
because of dynamical friction effects. As shown in Fig. [121, the position of r avoic j well corre- 
sponds to that of the observed minimum for all the clusters studied to date in a similar way. 
For NGC 6388, by adopting on = 18.9 km s" 1 and n n = 10 6 stars pc -3 as central velocity dis- 



persion and stellar density (jPryor fc Meylan Ill993f ), and by assuming the cluster structural 



parameters derived from the best-fit King model by lLanzoni et al. I (j2007d . see Table [2]), we 
have obtained r avoid ~ 15 r c . This is about 3 times larger than the location of observed min- 
imum, thus representing the first case where the two distances do not coincide. This result 
is quite puzzling and somehow suggests that NGC6388 appears "dynamically younger" than 
expected on the basis of its structural properties. In fact, our observations suggest that the 
dynamical friction in this cluster has been effective in segregating BSS (and similar massive 
objects) out to only 4-5 r c , whereas the theoretical expectation indicates that, within 15 r c , 
all stars with the mass of BSS or their binary progenitors should have already sunk to the 
center. Note that significantly larger (by a factor of 2) velocity dispersion, or lower (by a 
factor of 7) central density would be necessary to reconcile the expected and the observed 
minima. Why is dynamical friction less efficient in this cluster? Could the presence of an 
IMBH in the cluster center be important? As discussed in Lanzoni et al (2007), the radius 
of influence of a 5 x 10 3 M BH at the center of NGC 6388 is only 0.07 pc or 0.15 r c , so it 
is not obvious how it might affect cluster evolution at 5-15 r c . Perhaps the BSS that we 
observe at r > r c are stars which have "visited" the central region and have been put on 
highly eccentric orbits by the interaction with the BH. However this effect would probably 
fill-in the BSS avoidance region in the projected radial distribution. 

Could the BHk population, the suspected presence of a central IMBH, the anomalous 
position of the BSS radius of avoidance be connected to each other? While it is tempting to 
think that such a connection exists, we are at a loss to what that connection actually is. 
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Table 1. Papers in This Series 



NAME 



ARGUMENT 



REFERENCES 



M3 First detection of a BSS bimodal radial 

distribution 
M80 BSS population in the core 

47Tuc BSS radial distribution 

NGC288 BSS population in the core 
47Tuc BSS dynamical modeling and introduc- 

tion of zone of avoidance concept 
NGC6752 BSS radial distribution 
uj Cen The first and so far only flat BSS radial 

distribution 
47Tuc BSS chemical signatures 

NGC6266 BSS population in the core 

BSS dynamical modeling of four GCs 
M5 BSS obs and theoretical studies 

NGC 1904 BSS obs and theoretical studies 
M55 BSS radial distribution 



Ferraro et al. 1993,1997 

Ferraro et al. 1999 
Ferraro et al. 2003 
Bellazzini et al. 2002 
Mapelli et al. 2004 

Sabbi et al. 2004 
Ferraro et al. 2006a 

Ferraro et al. 2006b 
Beccari et al. 2006a 
Mapelli et al. 2006 
Lanzoni et al. 2007a 
Lanzoni et al. 2007b 
Lanzoni et al. 2007d 



Table 2. Adopted cluster parameters 



Center of Gravity a J20 oo = 17 h 36 m 17 s .23 <5j 20 oo = -44° 44' 07"1 

Concentration (c) 1.8 

Core Radius (r c ) 7"2 0.46 pc 

Tidal radius (r t ) 454" 29 pc 

Distance (m — M) = 15.60 13.2 kpc 

Reddening E(B - V) = 0.32 
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Table 3. The BSS population of NGC 6388 



Name 


RA 


DEC 


^255 


U 


B 


V 




[degree] 


[degree] 










BSS-1 


264.0713276 


-44.7367013 


19.586 


19.232 


19.995 




BSS-2 


264.0744196 


-44.7337510 


19.883 


19.252 


18.321 


17.876 


BSS-3 


264.0662077 


-44.7330923 


20.279 


19.453 


18.507 


17.842 


BSS-4 


264.0684376 


-44.7363620 


20.298 


19.463 


18.457 


17.709 



Note. — The complete table is available in electronic form. 



Table 4. Number Counts of BSS and HB sub-population Stars 





r e 


^BSS 


-Wrhb 


-^Vbhb 


^EHB 


A^BHk 


£samp I jjiaxap 


[arcsec] 


[arcsec] 


obs 


bck 


obs 


bck 


obs 


bck 


obs 


bck 


obs 


bck 







5 


22 





91 





18 





1 





6 





0.08 


5 


15 


43 





346 





57 





7 





8 





0.26 


15 


25 


17 





221 





43 





3 





5 





0.18 


25 


40 


10 


1 


247 


1 


54 





3 





4 





0.17 


40 


65 


24 


2 


248 


2 


52 


1 


8 





8 





0.16 


65 


110 


37 


7 


237 


6 


40 


1 


4 





6 





0.15 



Note. — The values listed as "obs" correspond to the observed number of stars counted, 
while those listed as "bck" quote the estimated number of contaminating field stars in each 
radial annulus (see Sect. 13.21) . 
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Fig. 1.— Map of the HST sample. The thick solid line delimits the WFPC2 HST FoV. The 
concentric annuli are used to study the radial distribution of BSSs. The inner and the outer 
annuli correspond to r = 5" and r = 110", respectively. 
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Fig. 2. — Map of the WFI sample. This sample has been used to estimate the structural 
parameters of NGC 6388 and the Galaxy contamination, but not for constructing the radial 
distribution of BSS. The dashed line marks the cluster tidal radius (r t = 454"). 
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Fig. 3. — Ultraviolet CMD of the HST WFPC2 sample. The different stellar populations 
discussed in the paper are marked with different symbols, as indicated by the labels. Solid 
dots mark the selected RR Lyrae stars. 
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Fig. 4.— Optical CMDs of the HST WFPC2 and the complementary ACS samples. 
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Fig. 5.— Comparison between the UV CMDs of NGC 6388 and M80. The latter has been 
suitably shifted in color and magnitude in order to superimpose the knees of the two HBs (at 
^255 — U ~ 0.2 and m 255 ~ 18.5). The dotted line marks the limit (771255 — 20) below which 
there are no more BHB stars in M80, and that we have adopted as the brightest boundary 
of the BHk population. 
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Fig. 6. — Zoom of the optical CMD of the complementary HST ACS sample, showing the 
different stellar populations discussed in the paper, marked as in Fig. [3j The selection box 
of BHk stars is marked as a dashed line, since this population is near the detection limits in 
the optical wavelengths. 
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Fig. 7. — Cumulative radial distributions of BHB, EHB and BHk stars (dashed lines), com- 
pared to that of the reference RHB population. No evidences of peculiar radial distributions 
are found for the blue HB sub-populations, with respect to normal cluster stars. 
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Fig. 8. — Differential Luminosity Function of NGC 6388 shifted to that of 47 Tuc. The 
dotted lines indicate the HB red clump and the RGB-bump level. From the inserts it is 
apparent that, other than the blue HB, the CMDs of the two clusters are quite similar. 
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Fig. 9. — Optical CMDs of the WFI sample for four different radial ranges, as marked by 
the labels in each panel. 
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Fig. 10. — Cumulative radial distribution of BSS (solid line) and RHB stars (dashed line) as 
a function of the projected distance from the cluster center, for the combined HST sample. 
The probability that they are extracted from the same population is ~ 10 -4 . 
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Fig. 11. — Radial distribution of the BSS and HB double normalized ratios, as defined in 
equation ( 14. 21) . plotted as a function of the radial coordinate expressed in units of the core 
radius. -Rrhb (with the size of the rectangles corresponding to the error bars computed 
as described in Sabbi et al. 2004) is almost constant around unity over the entire cluster 
extension, as expected for any normal, non-segregated cluster population. Instead, the radial 
trend of -Rbss (dots with error bars) is bimodal: highly peaked in the center (more than a 
factor of ~ 2 higher than Rrhb), decreasing at intermediate radii, and rising again outward. 
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Fig. 12. — Radial distributions of the specific fr equency Np,^/N m , a s observed for NGC 6388 
and for other four clusters (see references in lLanzoni et al.l 120071 ). The arrows mark the 
position of the estimated radius of avoidance (see Sect. [5]). This well corresponds to the 
position of the observed minimum of distributions, but in the case of NGC 6388. 



